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ABSTRACT 


In  this  report  we  complete  the  presentation  of  the  results  of  a  study  of  the  amplitudes 
of  the  quasi-static  electromagnetic  fields  produced  in,  on,  and  above  the  surface  of  a  sea 
of  finite  depth  by  vertically-directed  harmonic  electric  and  magnetic  dipoles  (VED’s  and 
VMD’s,  respectively)  and  horizontally-directed  harmonic  electric  and  magnetic  dipoles 
(HED’s  and  HMD's,  respectively).  In  a  previous  report  and  publication  we  presented 
numerical  data  for  the  fields  produced  on  and  above  the  surface  of  the  sea  by  »'ie  dipole 
sources  when  submerged  in  the  sea.  We  now  present  data  for  the  fields  produced  in  the 
sea  by  the  same  dipole  sources  when  they  are  either  on  or  elevated  above  the  sea  surface. 
The  primary  purpose  of  these  computations  is  to  determine  the  conditions  under  which 
an  electrically-conducting  sea  bed  can  produce  significant  changes  in  the  fields  generated 
by  the  dipoles,  as  compared  with  the  fields  produced  in,  on,  and  above  an  infinitely  deep 
sea,  for  frequencies  in  the  ULF,  ELF,  and  VLF  bands  (frequencies  less  than  30  kHz).  In 
contrast  to  our  earlier  work,  where  we  found  substantial  sea  floor  effects  on  the  fields 
produced  by  the  submerged  dipoles,  and  particularly  strong  effects  on  some  of  the  field 
components  produced  by  a  submerged  VED,  the  present  work  shows  that  the  fields 
produced  in  the  sea  by  elevated  dipoles  are  comparatively  little  affected  by  the  presence 
of  a  sea  floor.  This  result  is  in  fact  in  agreement  with  our  previously-derived  criteria  for 
estimating  when  a  sea  floor  effect  is  likely  to  occur:  (1)  the  sea  floor  effects  are  most 
marked  in  seas  that  are  electrically  shallow,  (2)  the  largest  changes  in  the  fields  occur 
when  the  dipoles  are  located  within  one  sea  water  skin  depth  of  the  sea  floor,  and  (3) 
there  are  essentially  no  sea  floor  effects  when  the  dipoles  are  located  more  than  three  sea 
water  skin  depths  above  the  floor.  These  new  data  and  results  should  prove  useful  in 
studies  of  air-undersea  communication  at  frequencies  less  than  30  kHz. 
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A  NOTE  ON  FREQUENCIES 


In  this  report  we  use  the  abbreviation  ULF  (ultra-low  frequencies)  for 
frequencies  less  than  5  Hz.  Pc  1  geomagnetic  pulsations  are  observed  in  the 
upper  part  of  this  frequency  range  (0.2  -  5  Hz).  ELF  (extremely-low 
frequencies)  is  used  to  designate  frequencies  in  the  range  5  Hz  to  3  kHz,  and 
VLF  (very-low  frequencies)  is  used  for  frequencies  in  the  range  3  to  30  kHz. 
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1.  INTRODUCTION 

This  report  extends  earlier  work  in  which  we  gave  a  broad  description  of  the  ultra-low 
frequency  (ULF;  frequencies  less  than  5  Hz)  and  extremely-low  frequency  (ELF; 
frequencies  in  the  range  5  Hz  to  3  kHz)  electromagnetic  fields  produced  on  and  above  the 
surface  of  a  sea  of  Finite  depth  by  submerged  harmonic  electric  and  magnetic  dipole 
sources  [Fraaer-Smith  and  Bubenik,  1979,  1080).  As  we  reported  at  the  time,  the 
purpose  of  this  earlier  work  was  to  provide  new,  quantitative  information  about  the 
effects  of  a  sea  floor  on  the  ULF /ELF  electromagnetic  fields  propagating  upwards  from 
submeiged  sources.  In  addition  to  providing  this  information,  our  work  showed  that  the 
sea  floor  effects  could  be  quite  large  under  certain  conditions,  particularly  at  frequencies 
in  the  ULF  range  and  for  a  vertical  electric  dipole  source.  Among  the  many  applications 
of  these  earlier  results,  we  pointed  out  two  specific  Navy  applications:  submarine 
detection  and  undersea-to-air  communication.  The  work  we  now  report  is  particularly 
relevant  to  the  latter  application;  our  results  describe  the  electromagnetic  field 
amplitudes  that  can  be  produced  in  the  sea,  and  the  effects  of  a  sea  floor  on  these 
amplitudes,  for  elevated  (i.e.,  airborne)  dipole  sources,  and  thus  the  data  are  applicable 
to  an  air-to-undersea  communication  link. 

When  combined  with  our  earlier  undersea-to-air  results,  these  new  air-to-undersea  data 
should  provide  much  of  the  quantitative  information  required  for  studies  of  two-way  air- 
undersea  communication  links  utilizing  frequencies  in  the  ULF,  ELF,  and  VLF  ranges. 
Since  radio  waves  with  frequencies  above  the  upper  limit  of  the  VLF  range  (30  kHz)  do 
not  penetrate  to  any  significant  depth  in  sea  water,  the  three  frequency  ranges 
specifically  covered  by  our  analysis  should  include  all  of  the  radio  frequencies  that  have 
any  practical  application  to  communication  with  moderately  or  deeply  submerged 
submarines. 

In  addition  to  their  communication  and  submarine  detection  applications,  our  earlier 
data  [Fraaer-Smith  and  Bubenik,  1979,  1980]  could  well  be  utilized  in  geophysical 
prospecting  and  in  the  derivation  of  effective  sea  floor  conductivities.  The  new  data 
presented  here  are  less  likely  to  have  application  in  these  latter  activities,  because  of  the 
much  smaller  relative  magnitudes  of  the  sea  floor  effects.  However,  there  is  a  possible 
application  analogous  to  submarine  detection,  namely  the  detection  by  electromagnetic 
means  of  the  presence  of  overhead  aircraft  by  submerged  submarines. 


Table  It  Representative  skin  depths  (for  t rg  =  4.0  S/m) 
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This  report  also  extends  earlier  work  in  which  we  considered  the  feasibility  of  air-to 
undersea  communication  at  ULF  using  airborne  loop  antennas  [Fraaer-Smith  et  at., 
1977].  The  frequencies  in  this  latter  study  were  largely  restricted  to  1  Hz  and  sea  floor 
effects  were  only  briefly  considered;  in  addition,  the  sources  of  the  ULF  fields  were 
lim.  id  to  loop  antennas,  to  magnetic  dipole  sources.  In  the  present  work  all  major 
categories  of  dipole  sources  are  covered  by  the  analysis:  vertically-directed  electric  and 
magnetic  dipoles  (VED’s  and  VMD’s,  respectively)  and  horizontally-directed  electric  and 
magnetic  dipoles  (HED’s  and  HMD’s,  respectively).  Further,  as  has  *>.lready  been 
mentioned,  all  frequencies  in  the  ULF/ELF/VLF  range  are  covered  and  information  is 
always  provided  about  the  possibility  of  sea  floor  effeeis.  Despite  the  limitations  of  its 
field  data,  the  earlier  work  is  not  superseded  since  its  primary  purpose  was  to  use 
representative  ULF  field  data  to  assess  the  feasibility  of  air-to-undersea  communication, 
whereas  the  purpose  of  the  present  work  is  to  present  a  wide  range  of  field  data  for  use 
in  such  studies. 

Our  computed  values  of  the  ULF/ELF/VLF  electromagnetic  fields  are  presented  in  a 
frequency-independent  parametric  form  that  has  two  features.  First,  all  distances  are 
measured  in  units  of  the  sea  water  skin  depth  S,  which  is  given  by 

6  —  (2/w/i0<r,)1/2,  (1.1) 

where  u  —  2,t /  is  the  angular  frequency,  m0  »s  the  permeability  of  free  space  (4*  X  10*' 
H/m),  and  a  is  the  electrical  conductivity  of  sea  water  measured  in  S/m.  The  second 
feature  is  the  use  of  a  sea  floor  conductivity  a j  that  is  measured  in  units  of  the  sea  water 
conductivity  <r  .  This  approach  enables  the  dipole  field  amplitudes  to  be  presented  in  a 
more  general  form  than  would  otherwise  be  possible  and  it  avoids  the  duplication  that 
would  result  from  the  use  of  specific  frequencies.  A  minor  disadvantage  of  the  approach 
is  the  need  for  a  computation  to  be  made  to  obtain  the  actual  electric  and  magnetic  field 
amplitudes  from  the  given  field  data,  which  are  also  in  parametric  form.  However,  these 
computations  are  comparatively  simple  and  can  be  made  quickly  on  a  hand  calculator; 
worked  examples  are  included  in  a  later  section  to  make  the  conversion  process  clear. 
The  same  parametric  approach  was  used  in  some  of  our  earlier  work  [Fraaer-Smith  and. 
Bubenik ,  1979,  1980],  and  further  worked  examples  are  given  there.  To  assist  further  in 
the  computation,  and  in  the  interpretation  of  the  data,  representative  skin  depths  for  sea 
water  (we  assume  rr}  =  4.0  S/m)  are  listed  in  Table  1. 

In  our  earlier  computations  for  the  undersea-to-air  path,  we  used  two  different  values 
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of  the  sea  floor  conductivity  to  show  the  effects  produced  by  tae  sea  floor  on  the  dipole 
fields  (specifically,  we  used  <r j  =  0.1<r4  and  <Tj  =  O.Cl^).  If  we  had  used  the  same 
approach  in  this  report,  we  would  have  had  to  include  nearly  100  figures,  which 
appeared  excessive.  We  therefore  restricted  our  computations  to  one  value  of  the  sea 
floor  conductivity,  namely  a j  =  0.1^.  While  this  restriction  certainly  limits  the 
quantity  of  our  data,  it  is  unlikely  to  seriously  limit  studies  using  the  data  because  of  the 
reduced  importance  of  sea  floor  effects  for  the  air-to-undersea  path.  This  is  consistent 
with  our  previous  observation  that  sea  floor  effects  are  most  pronounced  when  the  dipole 
sources  are  close  to  the  sea  floor  [Fraser-Smith  and  Bubenik ,  1979,  1980]. 

One  other  possible  limitation  of  our  data  results  from  the  restriction  of  the  dipole 
altitudes  to  the  three  values  h  —  0,  3.05,  and  10.05.  Because  our  work  is  primarily 
oriented  toward  air-to-undersea  communication  with  moderately  to  deeply  submerged 
receivers,  the  frequencies  of  most  interest  lie  in  the  ULF  and  lower-ELF  ranges,  where 
the  skin  depths  vary  from  a  few  tens  of  meters  up  to  many  kilometers  and  the  signals 
can  penetrate  to  considerable  depths  in  the  sea  water  before  they  become  too  attenuated 
to  be  detected.  For  these  skin  depths  the  restriction  of  A  to  a  maximum  value  of  105 
does  not  seriously  limit  the  range  of  applicability  of  our  data.  However,  for  frequencies 
in  the  VLF  range,  a  maximum  value  of  10.05  for  h  clearly  restricts  the  altitude  of  the 
dipole  sources  to  impractically  low  values  for  fixed-wing  aircraft  operation.  Fortunately, 
this  is  not  a  serious  limitation,  because  at  VLF  the  depth  of  any  sea  of  practical 
importance  for  air-to-undersea  communication  is  likely  to  exceed  10  -  100  skin  depths 
and  sea  floor  effects  can  be  ignored  (the  sea  water  skin  depth  varies  from  a  maximum  of 
about  4.6  m  at  3  kHz  down  to  1.4  m  at  30  kHz).  The  sea  can  then  be  assumed  to  be  of 
infinite  depth  and  the  VLF  fields  generated  beneath  its  surface  can  be  calculated  with 
little  difficulty  using  comparatively  simple  expressions  [e.g.,  Kraichman,  1976]. 

In  the  following  section  we  list  the  field  expressions  for  the  four  dipole  types  and 
describe  the  method  of  calculation  used  to  obtain  the  amplitudes  of  the  fields  (Section  2). 
We  next  briefly  describe  the  results  of  the  computations,  which  are  presented  almost 
entirely  in  graphical  form  in  six  Appendices  and  which  comprise  the  major  part  of  this 
report  (Section  3  and  Appendices).  Some  of  the  significance  of  these  results  are  discussed 
in  the  following  section  (Section  4).  A  list  of  references  (Section  5)  and  the  six 
appendices  complete  the  report.  In  general,  the  format  of  the  data  and  the  notation  are 
very  close  to  those  used  in  our  undersea-to-air  work  [  Fraser-Smith  and  Bubenik,  1979, 
1980].  However,  one  una>  oidable  difference  arises  because  of  the  interchange*  rce 


5 


and  Held  point  locations.  In  our  earlier  work  we  use  h  to  indicate  the  height  of  the  Held 
point,  or  receiver,  above  the  sea,  and  d  is  used  for  the  depth  of  the  dipole  source  in  the 
sea.  In  the  present  work  h  is  used  for  the  height  of  the  dipole  above  the  sea  and  d 
indicates  the  depth  of  the  field  point,  or  receiver,  in  the  sea.  All  the  other  symbols, 
including  D,  the  depth  of  the  sea,  have  the  same  meaning  as  before. 


Figure  Is  The  geometry  employed  in  computing  the  electromagnetic  field 
components  produced  at  point  P  located  at  a  depth  d  in  a  sea  of  finite  depth 
D  ( D>d )  by  dipole  sources  located  at  an  altitude  h  above  the  surface  of  the 
sea.  The  electromagnetic  field  components  for  two  dipole  sources  are 
illustrated  in  the  figure.  On  the  left  are  shown  the  Ef,  Eg  and  B ^  field 

components  produced  by  a  VED  source,  and  on  the  right  are  shown  the 
Bz,  B ^  and  Ef,  Ez;  E ^  components  produced  by  a  HMD  source. 
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2.  FIELD  CALCULATIONS 

Figure  1  shows  the  geometry  used  in  the  field  calculations;  except  for  the  interchange 
of  the  source  and  field  points,  it  closely  resembles  the  geometry  used  earlier  for  the 
undersea-to-air  link  [Fraser-Smith  and  Bubenik ,  1979,  1980).  A  cylindrical  coordinate 
system  (r,<t>,z)  is  used  and  the  harmonic  dipoles,  of  moment  m  (magnetic  dipole)  or  p 
(electric  dipole)  and  angular  frequency  u  (w  —  are  located  at  the  point  (0,0, |r\).  The 

sea  surface  is  the  plane  2  =  0,  the  sea  floor  is  the  plane  z  —  -D,  and  the  field 
measurements  are  made  at  the  point  P(r,^,-d),  where  D>d> 0.  The  region  2  >  0  is  air, 
assumed  to  be  free  space,  while  the  region  z  <  -£>  is  a  homogeneous  conducting  half¬ 
space  of  conductivity  oy  representing  the  sea  bed.  To  simplify  the  notation  in  the  field 
expressions,  the  air,  water,  and  sea  bed  regions  have  been  labeled  +1,  0,  and  -1, 
respectively,  and  the  surface  and  sea  floor  boundaries  are  correspondingly  denoted  by  +1 
and  -1. 

It  may  be  noted  that  the  HMD  source  in  Figure  1  is  directed  along  the  x  axis.  This 
particular  geometry,  which  is  also  used  for  the  HED,  differs  from  that  of  Kraxchman 
[1976],  who  oriented  his  HED  along  the  x-direction  and  his  HMD  along  the  y-direction. 

The  expressions  for  the  electromagnetic  fields  produced  by  the  VED,  VMD,  HED,  and 
HMD  at  the  point  P  were  derived  by  using  a  method  adapted  from  Morrison  et  ai, 
[1969],  and  Kong,  [1975].  The  same  method  was  used  in  our  earlier  work  to  obtain  the 
field  components  produced  above  the  sea  by  submerged  dipole  sources.  The  field 
expressions  for  the  four  dipole  types  are: 


1.  Vertical  Electric  Dipole 


Er  =  ~2 (2-1) 

(2>2) 

^  =  S  F™Xdy 


B, 


(2.3) 


2.  Vertical  Magnetic  Dipole 


(2.4) 

(2.5) 

(2.6) 

Er 

E* 

E, 


3.  Horizontal  Electric  Dipole 

p  cos  <b  T  Z00  1  /«= 

'—2^lk  'fo/awXdX-  -/o  {ioFiM  +  (\~io)F2E)<h 

psin  6  f  1  /*<•••■ 

"2  jrop“  [io  ^  x  dX  +  ^  /o  4M  +  (X  -  7n)  /^g)  A 

27T(fo  y<> 

POP  Sin  6  f  Z"00  70  „  1 

tT~  /«  ^-+>,F*E>'dx+-rl  f4£'),<x 


fl..= 


ft 


Pop  sin  <t> 


/: 


foE\'d\. 


2vr  •'0  70  +  X 

4.  Horizontal  Magnetic  Dipole 

£■  _  ^porrmin  <t>  f  1  Z00  1 

2/r  [r  Jo  70I 

= .4 


iwpom  cos  <j> 
2rr 


70  +  X 

1 


11.  .rFiE^dX 
|/°  7o  +  X  r 


X  <Aj, 

- 1  r 

r  Jo 


70 +  x 


Fog  X  (/X 


£  _  ^Wtwsin^ 

yc  F2MxrfX, 


(2.7) 

(2.8) 

(2.9) 

(2.10) 

(2.11) 

(2.12) 


(2-13) 

(2.14) 

(2.151 


Tiie  following  equations  also  apply 

F2M  =  F^Me-^Mh)M\r), 
FZM  =  F%e-('»d+>'h)M\r), 
FiM^F^Me-^d^M\r), 
F\e  -  Ffye-^^M^r), 
Foe  =  Ffte-^^^JUXr), 

FaE  =  F%c-^d+^M\r), 

FiE--=F{PEe-^d^Ji(\r), 

where 


Frlf  =  (I  +  P-\tm)IW  -  P-itmPitm) 
Fr\t  =  (1  -  P-ITAf)/U  -  P-lTAfPH'Ai) 
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^TE  =  (!  +  P-\Te)!{1  ~  P-ITEPITe)< 
FtE  =  (1  ~  P-lTE)/(l  ~  P-lTEPlTE )i 

pxtm  =  [m-n)/(r1  +  yo}]e"27orft 
P-iTM  =  [(Y-{  -  y0)/(r_,  +  %)]  t^-D\ 

PlTE  -  [(Zi  -  Z*)/{Zi  +  ft)]  e-2^, 

P-ire  =  [(*-1  -  Zo)/(^-i  +  Zo)\e2^d-D), 

Yn  —  (<rn  +  iuen)/%, 

Zn  =  iupKihn, 

*tn  -  X2  -  *a, 

*2  =  w2/i„«n  -  i<jjnn<rn. 


In  the  last  of  these  equations  the  subscript  n  can  take  one  of  the  three  values  1,  0,  or  -1, 
thus  specifying  the  applicable  region  for  the  subscripted  quantity.  Thus  <n  and  <rn,  for 
example,  are  the  permittivity  and  electrical  conductivity  in  region  n.  Note  that, 
according  to  this  notation  convention,  al  =  0,  <7Q  =  <rg,  and  —  oy.  The  equations  are 
valid  for  all  frequencies.  However,  in  our  computations  we  use  the  quasi-static 
approximation  (see  p.  3-15  of  Kraichman ,  (1076]),  which  is  equivalent  to  neglecting 
most  displacement  current  terms.  Thus,  with  the  single  exception  of  the  occurring  in 
the  expression  for  the  electric  field  component  Eg  produced  by  the  HMD,  we  set  «j  =  e0 
=  =  0.  This  is  a  commonly  used  approximation  that  is  applicable  when  the  source- 

receiver  distance  is  much  smaller  than  one  free  space  wavelength.  We  discussed  the 
frequency  range  over  which  this  approximation  is  valid  in  our  1980  report  [Fraser-Srmf/r 
and  Bubenik ,  1980]  and  concluded  that  the  field  data  only  begin  to  become  inaccurate 
for  frequencies  above  30  kHz.  The  same  arguments  used  in  that  earlier  report  are 
applicable  here  and  we  conclude  that  our  field  data  are  valid  throughout  the 
ULF/ELF/VLF  frequency  range. 

The  full  field  expressions  given  above  were  evaluated  numerically  using  the  same 
computational  techniques  used  in  our  previous  work  and,  in  particular,  using  the  same 
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techniques  described  by  Bubenik,  [1977].  Both  the  electric  dipole  moment  p  and  the 
magnetic  dipole  moment  m  are  set  equal  to  unity  (p  —  1  Am,  and  m  =  1  Am2).  For 
dipoles  of  arbitrary  moment  the  field  values  given  in  the  figures  should  be  multiplied  by 
the  moment  to  obtain  the  corresponding  field  magnitudes.  We  use  the  milUgamma  as 
our  unit  for  the  magnetic  field  (1  mq  =  1  picotesla),  and  the  electric  field  data  are 
presented  in  units  of  microvolts/meter. 


3.  RESULTS 


The  presentation  of  the  results  of  our  field  computations  follows  closely  the 
presentation  used  in  our  earlier  report  on  the  fields  generated  in  the  air  by  submerged 
dipoles  [Fraser-Smith  and  Bubenik,  lflSOj.  For  the  VED,  these  results  consist  of 
amplitude  data  for  the  radial  and  vertical  electric  field  components  E„  and  E.  the 
maximum  electric  field  vector  magnitude  EjQTAL,  and  the  maximum  magnetic  field 
vector  magnitude  5T0TaL  (which  represents  the  azimuthal  magnetic  field  component 
B ^).  For  the  VMD,  the  results  consist  of  amplitude  data  for  the  radial  and  vertical 
magnetic  field  components  Bf  and  Bg,  the  net  electric  field  component  ^qtal  (w^ch 
represents  the  azimuthal  electric  field  component  E^),  and  the  maximum  magnetic  field 
magnitude  5T0taL'  Also  presented  are  data  for  the  three  electric  field  components  (Ef, 
Ep  Eg),  three  magnetic  field  components  (Bf)  B ^  Bg),  the  maximum  magnetic  field 
magnitude  BT0Tal,  and  the  maximum  electric  field  magnitude  ^qtal  Pr°duced  by 
the  elevated  HED  and  HMD  along  the  two  principal  azimuthal  angles  ^  =  0°  and  00°. 
Although  it  is  not  the  primary  purpose  of  the  choice,  restriction  of  the  azimuthal  angles 
to  <f>  =  0°  and  00°  simplifies  the  presentation  of  field  data  for  the  two  horizontal  dipoles 
since  not  all  of  the  possible  field  components  are  produced  at  each  azimuthal  angle. 
Specifically,  only  one  of  the  horizontal  electric  components  (Er  and  E^)  or  one  of  the 
horizontal  magnetic  components  ( Br  and  B ^  are  produced  by  each  of  the  dipoles  when  <j> 
—  0°  or  90®,  and  the  vertical  field  components  ( E ,  and  B  \  also  vanish  at  one  or  the 

Z  Z' 

other  of  these  angles.  We  use  the  notation  ^HORIZONTAL  or  ^HORIZONTAL  ^or  non' 
zero  horizontal  field  components.  Whether  these  latter  components  are  radial  or 
azimuthal  can  be  determined  quickly  by  noting  the  azimuthal  angle  and  by  referring  to 
Equations  2,1  -  2.18.  Alternatively,  the  form  of  the  horizontal  components  can  be 
determined  by  reference  to  the  following  list  of  non-zero  electric  and  magnetic  field 
quantities  produced  by  each  dipole  type. 

1.  VED: 

•  The  basic  field  quantities  are  Er,  Eg,  B ^ 

•  Other  field  quantities  are  £T0TAL,  ^HORIZONTAL  *  Er’  *TOTAL  = 
^HORIZONTAL  =  Bp 


PREVIOUS  PAQE 
IS  BLANK 
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2.  VMD: 

•  The  basic  field  quantities  are  Bf,  Bg. 

•  Other  field  quantities  are  &TOTAL  =  ^HORIZONTAL  = 
^HORIZONTAL  “  B<t>'  ^TOTAL' 

3.  HED,  <t>  =  0°: 

•  The  basic  field  quantities  are  Efl  Eg,  B ^ 

•  Other  field  quantities  are  ^toTAL’  ^HORIZONTAL  =  Br>  BTOTAL  = 
^HORIZONTAL  “  Bf 

4.  HED,  4>  -  90°: 

•  The  basic  field  quantities  are  E ^  Br,  Bg. 

•  Other  field  quantities  are  fi^OTAL  *“  ^HORIZONTAL  =  ^TOTAL’ 
^HORIZONTAL  =  Br' 

5.  HMD,  *  =■•  0°: 

•  The  basic  field  quantities  are  E ^  Br,  Bg, 

•  Other  field  quantities  are  £T0TAL  =  ^HORIZONTAL  ~  B<t>'  ^TOTAL’ 
^HORIZONTAL  =  Br 

6.  HMD,  *  =  90°: 

•  The  basic  field  quantities  are  Er,  Ez,  B ^ 

•  Other  fie’d  quantities  are  ^tOTAL'  ^HORIZONTAL  “  Br'  "®TOTAL  = 
^HORIZONTAL  “  B<t>' 

Notice  that  there  are  only  three  basic  electric  and  magnetic  field  components  for  each 
dipole  category. 

The  presentation  of  field  data  is  very  similar  to  those  used  for  the  results  of  our 
submerged  dipole  computations  (Fraaer-Smi'f/j,  1079,  1980].  For  each  dipole  type  we 
first  present  a  series  of  curves  that  give  the  actual  fields  produced  by  the  elevated  dipole 
in  a  sea  of  infinite  depth,  i.e.,  without  any  effects  due  to  the  presence  of  a  sea  floor. 
Next,  we  present  additional  curves  that  give  the  ratios  of  the  fields  produced  by  the 
dipole  in  a  sea  of  finite  depth  D  to  the  fields  produced  under  otherwise  identical 
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conditions  in  the  sea  of  infinite  depth.  These  latter  curves  provide  an  immediate 
quantitative  indication  of  the  presence  or  absence  of  sea  floor  effects,  since  the  absence  of 
a  sea  floor  effect  is  indicated  by  a  ratio  of  1.0.  The  curves  for  the  infinitely  deep  sea  are 
presented  for  three  dipole  (t.e.,  transmitter)  altitudes  (h/S  =  0  (sea  surface),  3.0,  and 
10.0)  and  for  four  receiver  depths  (d/S  =  0.3,  1.0,  3.0,  and  10.0).  These  data  are 
supplemented  by  the  sea  floor  effect  curves,  which  are  drawn  for  four  sea  depths  (D/6  — 
0.3,  1.0,  3.0,  and  10.0).  Only  one  sea  bed  electrical  conductivity  is  considered;  the  value 
chosen,  c .  =  0.1  <r is  intended  to  be  typical  of  the  values  encountered  in  the  first  1  km 
depth  of  the  sea  bed.  The  reader  is  referred  to  the  work  of  Young  and  Cox,  [1081],  for 
the  results  of  some  recent  measurements  of  sea  bed  conductivity  that  support  our  choice. 
Note  that  the  subscript  has  been  dropped  from  a.  when  this  conductivity  appears  in  the 
ordinate  labels  of  some  of  the  infinitely-deep-sea  figures. 

The  field  data  are  presented  in  the  following  figures: 

1.  VED;  Figures  2  through  9  (Appendix  1). 

2.  VTVtD;  Figures  10  through  17  (Appendix  2). 

3.  HED,  <p  =s  0°;  Figures  18  through  25  (Appendix  3). 

4.  HED,  <t>  —  90°;  Figures  26  through  33  (Appendix  4). 

5.  HMD,  <6  =  0°;  Figures  34  through  41  (Appendix  5). 

6.  HMD,  i  aw.  90°;  Figures  42  through  49  (Appendix  6). 

The  same  basic  arrangement  of  figures  is  used  for  each  dipole  category.  One  small 
difference  occurs  in  the  case  of  the  HMD  (^  =  90°),  where  fewer  sea-floor-effect  curves 
are  shown.  This  difference  is  caused  by  the  comparatively  small  amplitude  of  the  E , 
component,  Since  ^HORIZONTAL  =  ^TOTAL’  a  c^ose  approximation,  the 
^HORIZONTAL  curves  have  been  omitted, 

To  illustrate  the  use  of  these  data,  suppose  the  source  of  fields  is  an  HED  ( <t>  =  90°) 
transmitting  at  1  Hz  and  we  wish  to  know  Bs  at  a  distance  r/6  =  100  and  depth  d/i  = 
0.3  for  a  dipole  altitude  h/6  =  10.  Our  B  curves  for  the  HED  (<£  =  90°)  give  B  X  6  = 
3.8  X  10'2  m7.m2  (Figure  27).  At  1  Hz,  6  —  251.6  in  and  Br  —  3.0  X  10'7  m7  for  a  unit 
moment  dipole  and  infinitely  deep  sea.  Turning  next  to  Figure  33,  we  find  that  the  sea 
floor  effect  (for  cj  —  0.4  S/m)  varies  from  negligible  (multiplication  factor  ■=  1)  for  a  sea 
depth  of  1.0  or  more  sea  water  skin  depths,  i.e.,  for  depths  greater  than  about  250  m, 
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4.  DISCUSSION 

In  our  studies  of  the  ULF/ELF  electromagnetic  fields  generated  above  a  sea  of  finite 
depth  by  submerged  harmonic  dipoles  [Fraaer-Smith  and  Bubenik,  1979,  1980]  we  drew 
attention  to  the  differing  power-law  declines  with  distance  of  the  field  components 
produced  by  the  various  dipole  sources  at  large  distances.  The  purpose  of  this 
comparison  was  to  identify  the  dipole  types  and  field  components  with  the  greatest 
potential  for  undersea-to-air  communication.  We  finally  concluded  that  the  E ,  field 
produced  by  an  HED  (<j>  =  0°)  was  likely  to  be  the  most  attractive  electric  field 
component.  Several  different  factors  led  to  this  conclusion,  among  which  the  inverse 
square  decline  of  the  amplitude  of  the  field  component  with  distance  (at  large  distances) 
was  perhaps  the  most  influential.  We  now  examine  the  large-distance  variations  of  the 
various  field  components  considered  in  the  present  study,  together  with  the  effects  of  a 
sea  floor  on  these  components,  with  the  goal,  once  again,  of  determining  which  dipole 
types  and  field  components  are  likely  to  have  the  greatest  potential  for  communication 
over  long  distances. 

At  horizontal  distances  greater  than  105  the  curves  for  the  ULF/ELF  fields  produced  in 
an  infinitely  deep  sea  by  our  elevated  dipole  sources  begin  to  decrease  linearly  with 
distance  (see  Figures  2,  3  [VED];  10,  11  [VMD];  18,  19  [HED,  =  0°];  26,  27  [HED,  <t>  = 
90°];  34,  35  [HMD,  <t>  —  0°];  42,  43  [HMD,  <t>  =  90°]).  Because  of  the  use  that  is  made  of 
logarithmic  scales,  the  linear  decline  implies  a  power  law  decrease  of  the  fields  with 
distance:  the  fields  decrease  with  distance  as  r‘n,  where  the  power  n  depends  both  on 
the  dipole  type  and  on  the  electric  or  magnetic  field  component  under  consideration.  As 
indicated  above,  power  law  decreases  of  this  kind  are  also  observed  in  the  fields 
produced  in  the  air  by  submerged  dipoles,  Further,  they  are  in  agreement  with  the 
power  law  variations  in  approximate  expressions  for  the  fields  that  apply  at  large  source- 
receiver  distances  [ Kraichman ,  1076;  also  see  Bannister,  1984,  for  a  new  series  of 
formulas  that  apply  at  measurement  distances  greater  than  105  from  the  source]. 

The  values  of  n  for  the  various  field  components  produced  in  an  infinitely  deep  sea  by 
our  various  categories  of  elevated  dipoles  are  listed  in  Table  2.  It  can  be  seen  that  n 
ranges  from  5,  implying  a  very  rapid  decline  of  the  field  quantity  with  distance,  down  to 
2,  implying  a  comparatively  slow  variation  of  the  field  quantity  with  distance.  The 
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Table  2t  Values  of  the  power  n  occurring 
r'n  observed  in  the  fields  produced  in  the  sea 
asterisk  indicates  that  no  field  component  is  ] 


VMD 

VED 

HED, 

4  =  o9 

5T0TAL 

4 

2 

3 

^TOTAL 

4 

2 

3 

Br 

4 

* 

* 

B* 

* 

2 

3 

Bz 

5 

* 

* 

Er 

* 

2 

3 

E, * 

4 

* 

* 

Ez 

* 

3 

4 
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distance  variations  corresponding  to  these  powers  of  n  are  all  in  agreement  with  the 
variations  predicted  by  the  approximate  expressions  for  the  fields  tabulated  by 
Kraichman  [1976]  for  large  source-receiver  distances. 

Our  primary  interest  is  in  long-distance  air-to-undersea  communication,  for  which  field 
quantities  having  a  r 2  distance  variation  are  likely  to  be  most  suitable.  However,  we 
recognize  that  there  may  be  situations  where  a  short  communication  range  could  be  an 
important  requirement.  Such  a  situation  might  arise,  for  example,  during  a  naval 
engagement,  where  good  communication  would  be  important  on  a  local  basis  but  where 
it  would  also  be  particularly  desirable  that  the  enemy  fleet  (presumably  more  distant) 
should  be  denied  access  to  the  communication  signals.  For  such  limited-range 
communication  the  VMD  appears  to  be  the  most  suitable  transmitter,  since  none  the 
fields  it  produces  have  an  n  less  than  4  (Table  2),  whereas  all  the  other  dipole  sources 
produce  fields  in  the  sea  with  a  smaller  n. 

Examining  the  field  quantities  with  n  =  2  in  greater  detail,  we  observe  that  they  are 
Er  (or,  equivalently,  ^horizontal)  and  D +  (or'  equivalently,  flT0TAL)  for  the  VED, 
and  E2  for  the  HMD  (<j>  —  90°).  These  three  field  components  would  appear  to  have 
equivalent  potential  for  long  distance  communication,  but  we  can  conditionally  dismiss 
the  Ex  component  produced  by  the  HMD  source  because  its  amplitude  is  comparatively 
small  (over  the  distance  range  considered  in  our  computations  it  is  smaller  than  the  E ^ 
component  produced  by  the  same  source,  even  though  the  amplitude  of  the  latter 
component  declines  more  rapidly  with  distance).  This  leaves  the  VED  as  the  preferred 
source  for  long-distance  ai» -to- undersea  communication. 

The  conditionality  of  our  dismissal  of  the  HMD  source  for  long-distance  communication 
follows  from  the  difficulty  of  comparing  the  fields  produced  by  electric  and  magnetic 
dipole  sources.  We  assume  unit  moment  dipolo  sources  wheix  comparing  the  magnitudes 
of  the  field  quantities.  However,  there  are  many  practical  differences  between  the 
sources,  which  influence  the  magnitude  oT  the  dipole  moment  that  can  be  achieved,  and 
technological  advances  can  also  alter  the  situation.  At  present  there  is  a  possibility  that 
a  very  powerful  HMD  source  might  be  developed  (i.e.,  a  rotating  superconducting  loop). 
If  such  a  development  were  to  occur,  our  conclusion  that  the  VED  is  the  preferred  source 
for  long-distance  communication  would  need  to  be  reexamined. 

In  our  earlier  work  [Eraser-Smith  and  Bubenik,  1979,  1980]  we  found  that  the  fields 
produced  above  the  sea  by  a  VED  were  subject  to  exceptionally  large  sea  floor  effects: 


under  certain  conditions  the  Er  field  component  can  be  increased  by  nearly  six  orders  of 
magnitude,  Such  large  effects  are  not  observed  in  our  results  for  the  elevated  VED.  Our 
computations  indicate  that  field  enhancements  by  a  factor  of  two  are  about  the  best  that 
can  be  expected  at  larger  distances  {&roTAL  and  ^HORIZONTAL,)  an<*  Auctions  °f  fold 
amplitudes  to  about  half  their  original  amplitude  can  also  take  place  (^-poTAL  an<* 

The  sea  floor  effects  are  therefore  not  particularly  marked  and  do  not  influence  our 
conclusion  regarding  the  advantages  of  a  VED  as  a  transmitter. 

Taken  as  a  whole,  our  computations  show  that  the  fields  produced  by  ail  the  dipole 
sources  are  little  affected  by  the  presence  of  a  sea  floor  except  for  the  Et  component 
produced  by  the  HED  {<j>  —  08).  In  this  latter  case  very  substantial  enhancements  can 
occur  for  horizontal  distances  in  the  range  5  <  r/6  <  100,  provided  the  dipole  is  located 
on  the  sea  surface.  However,  once  the  dipole  is  elevated  above  the  sea  surface  the 
enhancements  rapidly  decline  to  the  same  order  as  those  for  the  other  dipole  categories. 
It  does  not  appear  that  sea  floor  effects  are  likely  to  be  a  substantial  factor  in  air-to- 
undersea  communication. 

It  is  interesting  to  observe  that  the  criteria  we  developed  in  our  submerged  dipole 
studies  [Fraaer-Smith  and  Bubenik,  1070,  1080]  for  estimating  when  sea-floor  effects  are 
most  likely  to  occur  are  also  valid  for  the  ULF/ELF/VLF  fields  generated  in  the  sea  by 
elevated  dipoles:  (1)  the  sea-floor  effects  are  most  marked  in  seas  that  are  electrically 
shallow,  (2)  the  largest  changes  in  the  fields  occur  when  the  dipoles  are  located  within 
onesea  water  skin  depth  of  the  sea  floor,  and  (3}  there  are  essentially  no  sea-floor  effects 
when  the  dipoles  are  located  more  than  threesea  water  skin  depths  above  the  floor. 
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Field  Data  For  The 
VERTICAL  ELECTRIC  DIPOLE 
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tically  directed  harmonic  electric  dipole  (VED)  located  on  and  above  the  surface  of 
the  sea. 


Figure  3.  Variation  with  horizontal  distance  of  the  amplitudes  of  the  two  electric 
field  components  ^HORIZONTAL  an<*  Ez  produced  in  an  infinitely  deep  sea  by  a  verti¬ 
cally  directed  harmonic  electric  dipole  IVED)  located  on  and  above  the  surface  of 
the  sea. 


electrically  conducting  sea  floor  (conductivity  =  0.1  os  ,  depth  D)  is 
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Figure  7.  Curves  illustrating  the 
data  presented  in  Figure  3  for  a 
surface  of  an  infinitely  deep  sei 
floor  (conductivity  of  =  0.1  as> 


(BTOTAL)  and  electric  (EyoTAL)  fields  produced  in  an  infinitely  deep  sea  by  a  veri 
cally  directed  harmonic  magnetic  dipole  (VMD)  located  on  and  above  the  surface  of 
the  sea. 


data  presented  in  Figure  10  for  a  VHD  located  or.  the  surface  of  an  inf'.mtsiy  deep 
sea  when  an  electrically  conducting  sea  floor  (conductivity  of  -  0.1  os,  depth  D)  is 
introduced. 
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when  an  electrically  conducting  sea  floor  (conductivity  Of  =  0.1  os,  aepth  D) 
n traduced. 
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Field  Data  For  The 

HORIZONTAL  ELECTRIC  DIPOLE,  0° 
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(Bjotal)  and  electric  (EtotAl)  fields  produced  in  an  infinitely  deep  sea  fay  a  hori¬ 
zontally  directed  harmonic  electric  dipole  (HED)  located  on  and  above  the  surface  of 
the  sea.  The  fields  are  given  for  an  azimuthal  (<(>)  of  0°. 
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21.  Curves  illustrating  the  changes  produced  in  the  electric  field  component 
presented  in  Figure  19  for  an  HED  located  on  the  surface  of  an  infinitely 
sea  when  an  electrically  conducting  sea  floor  (conductivity  of  =  0.1  as,  depth 
introduced.  Note  that  the  changes  apply  to  the  fields  produced  at  an  azimuthal 


above  the  surface  of  an  infinitely  aeep  sea  vn/6  =  wnen  an  eiectricaiiy 
conducting  sea  floor  (conductivity  of  =  0.1  o5,  depth  D)  is  introduced.  Note 
that  the  changes  apply  to  the  fields  produced  at  an  azimuthal  angle  (<f)  of  0°. 
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Figure  23.  Curves  illustrating  the  changes  produced  in  the  electric  field  component 
data  presented  in  Figure  19  for  an  HED  located  three  sea  water  skin  depths  above 
the  surface  of  an  infinitely  deep  sea  (h/6  =  3.0)  when  an  electrically  conducting 
sea  floor  (conductivity  af  =  0.1  05,  depth  D)  is  introduced.  Note  that  the  changes 
apply  to  the  fields  produced  at  an  azimuthal  angle  (<J>)  of  0°. 
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Field  Data  For  The 

HORIZONTAL  ELECTRIC  DIPOLE,  </»—  00° 
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zontally  directed  hanucnic  electric  dipcle  (HED)  located  on  and  above  the  surface  of 
the  sea.  The  fields  are  given  for  ar.  eziauthal  angle  («J>)  of  90°. 


data  presented  in  Figure  26  for  an  HED  located  on  the  surface  of  an  infinitely  deep 
sea  when  an  electrically  conducting  sea  floor  (conductivity  Of  =  0.1  os,  depth  D)  is 
introduced.  Note  that  the  changes  apply  to  the  fields  produced  at  an  azinuthal  angle 


Figure  31 .  Curves  illustrating  the  changes  produced  in  the  magnetic  field  component 
data  presented  in  Figure  27  for  an  HED  located  three  sea  water  skin  depths  above 
the  surface  of  an  infinitely  deep  sea  (h/6  =  3.0)  when  an  electrically  conducting 
sea  floor  (conductivity  op  =  0.1  crs,  depth  D)  is  introduced.  Note  that  the  changes 
apply  to  the  fields  produced  at  an  azimuthal  angle  ($)  of  90°. 


field  data  presented  in  Figure  26  for  an  HED  located  ten  sea  water  skin  depths 
above  the  surface  of  an  infinitely  deep  sea  (h/6  =  10.0)  when  an  electrically 
conducting  sea  floor  (conductivity  Of  =  0.1  as,  depth  D)  is  introduced.  Note 
that  the  changes  apply  to  the  fields  produced  at  an  azimuthal  angle  ($)  of  90°. 
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gure  33.  Curves  illustrating  the  changes  produced  in  the  magnetic  field  component 
data  presented  in  Figure  27  for  an  HED  located  ten  sea  water  skin  depths  above  the 
surface  of  an  infinitely  deep  sea  (h/6  =  10. 0)  when  an  electrically  conducting  sea 
floor  (conductivity  of  =  0.1  crs,  depth  D)  is  introduced.  Note  that  the  changes 
apply  tG  the  fields  produced  at  an  azinuthal  angle  (4?)  of  90°. 
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Field  Data  For  The 


HORIZONTAL  MAGNETIC  DIPOLE,  0  ~  0° 
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zontally  directed  harmonic  aagnetic  dipole  (HMD)  located  on  and  above  the  surface  of 
the  sea  The  fields  are  given  for  an  aziRuthal  angle  (d>)  of  0°. 
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Figure  35.  Variation  with  horizontal  distance  of  the  amplitudes  of  the  two  magnetic 
field  components  BhqrizoNTAL  and  Bz  produced  in  an  infinitely  deep  sea  by  a  hori¬ 
zontally  directed  harmonic  magnetic  dipole  (HMD)  located  on  and  above  the  surface 
of  the  sea.  The  fields  are  given  for  an  azimuthal  angle  (<J>)  of  0°. 
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surface  of  an  infinitely  deep  sea  (h/6  =  3.0)  when  an  electrically 
a  floor  (conductivity  af  =  0.1  crs,  depth  D)  is  introduced.  Note  t! 
■s  apply  to  the  fields  produced  at  an  azimuthal  angle  (<f>)  of  0°. 


data  presented  in  Figure  35  for  an  HMD  located  three  sea  water  skin  depths  above 
the  surface  of  an  infinitely  deep  sea  (h/6  =  3.0}  when  ar.  electrically  conducting 
sea  floor  (conductivity  of  =  0.1  os,  depth  C)  is  introduced.  Hate  that  the  changes 
apply  to  the  fields  produced  at  an  azimuthal  angle  {<{>)  of  0°. 


above  the  surface  of  an  infinitely  deep  sea  (h/5  =  10.0)  when  an  electrically 
conducting  sea  floor  (conductivity  Of  =  0.1  as,  depth  D)  is  introduced.  Note 
that  the  changes  apply  to  the  fields  produced  at  an  azimuthal  angle  ($)  of  0° 
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APPENDIX  5 


Field  Data  For  The 

HORIZONTAL  MAGNETIC  DIPOLE,  <£*=  00° 


(BtotAL)  and  electric  (EjOTAL)  fields  produced  in  an  infinitely  deep  sea  by  a  hori¬ 
zontally  directed  harmonic  magnetic  dipole  (HMD)  located  on  and  above  the  surface  of 
the  sea.  The  fields  are  given  for  an  azimuthal  angle  (a)  of  90°. 
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field  data  presented  in  Figure  42  for  an  HMD  located  on  the  surface  of  an  infi¬ 
nitely  deep  sea  when  an  electrically  conducting  sea  floor  (conductivity  Of  =  0. 
depth  D)  is  introduced.  Note  that  the  changes  apply  to  the  fields  produced  at  , 
azimuthal  angle  (g>)  of  90°. 


DISTRIBUTION  LIST 


No.  of 

Organization  Copies 


No.  of 

Organization  Copies 


Director 

Defense  Advanced  Research 
Projects  Agency 

ATTN:  Project  Management  2 

GSD,  R.  Alewine  1 

STO,  D.D.  Lewis  1 

1400  Wilson  Blvd 
Arlington,  VA  22200 


Defense  Technical 

Information  Center  12 

Cameron  Station 
Alexandria,  VA  22314 

Office  of  Naval  Research 
ATTN:  Code  222  1 

Code  414  1 

Code  420  1 

Code  425  1 

800  North  Quincy  Street 
Arlington,  VA  22217 


Office  of  Naval  Research 

Resident  Representative  1 

University  of  California, 

San  Diego 
La  Jolla,  CA  02093 

Office  of  Naval  Research 
Resident  Representative  1 

Stanford  University 
Durand  Building,  Room  165 
Stanford,  CA  94305 

Assistant  Deputy  Undersecretary 
of  Defense  (C3) 

ATTN:  T.P.  Quinn  1 

Pentagon 

Washington,  D.C  20301 


Naval  Research  Laboratory 
Information  Technology  Division 


ATTN:  J.R.  Davis  1 

4555  Overlook  Avenue,  3W 
Washington,  D.C.  20375 

Naval  Ocean  Systems  Center 
ATTN:  Library  1 

K.L.  Grauer  1 

C.F.  Ramstedt  1 

Y.  Richter  1 

271  Catalina  Boulevard 
San  Diego,  CA  95152 


Naval  Electronic  Systems 
Command 

ATTN:  PME-110-112  1 

PME-110-X1  1 

Department  of  the  Navy 
Washington,  D.C,  20360 

Naval  Underwater  Systems  Center 


New  London  Laboratory 
ATTN:  P.  Bannister  1 

A.  Bruno  1 

J.  Orr  1 

E.  Soderberg  1 

New  London,  CT  06320 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
ATTN:  J.J.  Holmes  1 

P.  Wessel  I 

Silver  Spring,  MD  20910 

David  W.  Taylor  Naval  Ship 
Research  &  Development  Center 
ATTN:  W.  Andahazy  1 

P.  Field  1 

Annapolis,  MD  21402 


DISTRIBUTION  LIST  (Continued) 


No.  of 

No.  of 

Organization 

Copies 

Organization 

Copies 

Office  of  the  .Assistant  Secretary 

Naval  Air  Development  Center 

of  the  Navy  (R,  E&S) 

ATTN:  J.  Shannon 

1 

ATTN:  J.  Hull 

1 

Warminster,  PA  18974 

Washington,  D.C.  20350 

Naval  Ocean  Systems  Center 

Naval  Ocean  R  &  D  Activity 

ATTN:  R.  Buntzen 

1 

ATTN:  D.L.  Durham 

1 

C.  Fuzak 

1 

D.W.  Handschumacher 

1 

P.  Hansen 

1 

K.  Smits 

1 

271  Catalina  Boulevard 

NSTL  Station 

Bay  St.  Louis,  MS  39522 

San  Diego,  CA  95152 

Director 

Naval  Oceanographic  Office 

Defense  Nuclear  Agency 

ATTN:  O.E.  Avery 

1 

ATTN:  RAAE 

1 

T.  Davis 

1 

DOST 

1 

G.R.  Lorentzen 

1 

RAEV 

1 

NSTL  Station 

Washington,  D.C.  20305 

Bay  St.  Louis,  MS  39522 

SRI  International 

Naval  Air  Systems  Command 

ATTN:  D.M.  Bubenik 

1 

ATTN:  B.L.  Dillon 

1 

J.B.  Chown 

1 

Washington,  D.C.  20361 

J.G.  Depp 

1 

R.C.  Honey 

1 

Naval  Intelligence 

W.E.  Jaye 

1 

Support  Center 

333  Ravenswood  Avenue 

ATTN:  G.D.  Batts 

1 

Menlo  Park,  CA  94025 

W.  Reese 

1 

4301  Suit  land  Road 

R.  &  D.  Associates 

Washington,  D.C.  20390 

ATTN:  C.  Greifinger 

P.O.  Box  9695 

1 

Naval  Postgraduate  School 
Department  of  Physics 

Marina  del  Rey,  CA  90291 

and  Chemistry 

Pacific-Sierra  Research  Corp. 

ATTN:  O.  Heinz 

1 

ATTN:  E.C.  Field,  Jr. 

1 

Monterey,  CA  93940 

12340  Santa  Monica  Blvd 

Los  Angeles,  CA  90025 

Naval  Coastal  Systems  Center 
ATTN:  R.H.  Clark 

TASC 

M.J.  Wynn 

1 

ATTN:  J.  Czika,  Jr. 

1 

Panama  City,  FL  32407 

1700  N.  Moore  Street,  Suite  1220 
Arlington,  VA  22209 

ii 


REPRODUCED  FROM 
BEST  AVAILABLE  COPY 


DISTRIBUTION  LIST  (Continued) 


No.  of 

Organization  Copies 


Naval  Weapons  Center 

ATTN:  R.J.  Dinger  1 

China  Lake,  CA  93555 

Johns  Hopkins  University 
Applied  Physics  Laboratory 
ATTN:  L.W.  Hart  1 

H.  Ko  1 

Johns  Hopkins  Road 
Laurel,  MD  20810 

La  Jolla  Institute 

ATTN:  K.  Watson  1 

La  Jolla,  CA  92407 

E.G.  &  G. 

ATTN:  L.E.  Pitts  1 

P.O.  Box  398 
Riverdale,  MD  20840 


University  of  Texas,  Austin 
Ge^magnetics  and  Electrical 
Geoscience  Laboratory 
ATTN:  F.X.  Bostick,  Jr.  1 

Austin,  TX  78712 

Commander 
Rome  Air  Development 
Center/EEPS 

ATTN:  P.A.  Kossey  1 

Hanscom  AFB,  MA  01731 


No:  of 

Organization  Copies 


University  of  California 
Scripps  Institute  of  Oceanography 
ATTN:  C.S.  Cox  1 

La  Jolla,  CA  92093 

Lockheed  Palo  Alto  Research 
Laboratory 

ATTN:  J.B.  Cladis  1 

W.L.  Imhof  1 

J.B.  Reagan  1 

M.  Walt  1 

3251  Hanover  Street,  Bldg  255 
Palo  Alto,  CA  94304 


Chief 

Air  Force  Technical 
Applications  Center  1 

HQUSAF 

Patrick  AFB,  FL  32925 
Commander 

Air  Force  Systems  Command  2 

Andrews  AFB,  MD  20331 

Commander 

Air  Force  Geophysics  Laboratory 
ATTN:  Technical  Library  1 

Hanscom  AFB,  MA  01731 


REPRODUCED  FROM 
BEST  AVAILABLE  COPY 


